Doped Spin Liquid: Luttinger Sum Rule and Low Temperature Order 
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We analyze a model of two-leg Hubbard ladders weakly coupled by interladder tunneling. At 
half filling a semimetallic state with small Fermi pockets is induced beyond a threshold tunneling 
strength. The sign changes in the single electron Green's function relevant for the Luttinger Sum 
Rule now take place at surfaces with both zeroes and infinities with important consequences for the 
interpretation of ARPES experiments. Residual interactions between electron and hole-like quasi- 
particles cause a transition to long range order at low temperatures. The theory can be extended 
to small doping leading to superconducting order. 



While the properties of doped spin liquids in more than 
one dimension are notoriously difficult to analyze, they 
are nonetheless highly relevant. In one dimension the two 
leg Hubbard ladder at half filling is the spin liquid epit- 
ome and as such ladder systems have attracted strong 
interest (for an early review see 0). Powerful analytic 
techniques such as bosonization and Bethe ansatz have 
been applied to single ladders with weak interactions (see 
and references therein) and have led to comprehen- 
sive understanding of both doped and undoped ladders. 
In this letter we report the extension of these results to 
higher dimensions through the introduction of a small 
long range interladder tunneling in an ensemble of un- 
coupled half-filled Hubbard ladders. Increasing the tun- 
neling amplitude leads to the formation of closed elec- 
tron and hole Fermi pockets. The Luttinger Sum Rule 
(LSR) now takes on a novel form with the sign changes 
in the one electron Green's function appearing both as 
zeroes and infinities. This result has strong implications 
for the interpretation of ARPES results on underdoped 
cuprates. In the pseudogap phase the experiments in- 
terpret infinities as a set of disconnected Fermi arcs 0, 
but do not (and cannot) observe the zeroes. Lastly we 
analyze possible instabilities of the carriers in the Fermi 
pockets using interactions derived from the low energy 
effective field theory for the ladders. 

The dynamics of the component half-filled ladders in 
our ensemble are governed at low energies by an effective 
field theory. As demonstrated in 0, half-filled ladders 
with gencrically repulsive interactions experience under 
renormalization an enhancement in the symmetry of the 
bare Hubbard lattice Hamiltonian. With this enhance- 
ment, the effective field theory for the ladder is the SO (8) 
Gross- Neveu model, H so ^ [4J. We couple the half-filled 
ladders together by long range single-particle tunneling. 
The complete Hamiltonian is then 



SO(8) 



x)c n i'i' (T (x) + h.c.),(l) 



amplitude long range, we acquire a small parameter, as 
was done for a similar model of coupled Hubbard chains 
in . We assume the following hierarchy of energy scales: 
W (bandwidth) > A (spectral gap) > t 1 - . The first 
inequality guarantees that the ladders can be described 
using the effective field theory. 

The Gross-Neveu model is exactly solvable for all semi- 
simple symmetry groups and a great deal is known about 
its thermodynamics and correlation functions. In the 
SO (8) case the correlation functions were studied in 0,0]. 
The spectrum of this model consists of three octets of 
particles of mass A and a multiplet of 29 excitons with 
mass y/3 A. Two octets consist of quasi-particles of dif- 
ferent chirality transforming according to the two irre- 
ducible spinor representations of SO (8), while the third 
octet consists of vector particles. The latter include mag- 
netic excitations as well as the Cooperon (a particle with 
charge ±2e). The 16 kink fields, carrying charge, spin, or- 
bit, and parity indices, are direct descendants of the orig- 
inal electron lattice operators on the ladders. According 
to , the corresponding single electron Green's func- 
tion is given by 



z„ 



P - el(k x ) - A^ 



(2) 



where c^ H(T creates an electron at the n-th site on the i-th 
leg (I = 1, 2) of the i-th ladder. By making the hopping 



where a = a,j with a spin and j = ± indexing the 
bonding/antibonding bands. e a (k x ) is the bare disper- 
sion in the corresponding band. There are no off diagonal 
Green's function involving electrons from opposite Fermi 
points, a result of right and left moving quasi-particles 
belonging to different irreducible representations of the 
SO(8) group. As was demonstrated in 0, the incoherent 
part of the Green's function G a ,reg yields a negligible con- 
tribution to the spectral weight. Thus the quasi-particle 
weight satisfies Z„«l. 

RPA Analysis: We will study the properties of our cou- 
pled ladders close to the Mott-Hubbard transition. Our 
approach follows closely the one developed in . The in- 
terladder hopping is diagonal in the bonding/antibonding 
indices because of the bands' differing Fermi wavevectors. 
A Random Phase approximation (RPA) (diagrammati- 
cally pictured in Fig. 1) yields the following expression 
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FIG. 1: The RPA equation for single particle Green's func- 
tion (thick line), Ga PA ■ The double line is the bare hopping 
amplitude, t a while the thin line is the bare Green's function 
G . 



for the full 2D Green's function: 



G 



RPA 



(fer,kj_) = 



{{G^iK))- 1 -^)}' 1 . 



(3) 



In our model the interchain tunneling amplitude has 
strong peaks at kj^ = 0, G/2, where G is the inverse lat- 
tice vector in the direction perpendicular to the chains. 
(The peak at G/2 follows from particle- hole symmetry, 
i.e. t^(k) = + G/2).) Near these points the fol- 

lowing expansion is valid: 



*„(kx + (IT l)G/4) = T4 a o[l - (k±) 2 /^o 



(4) 



where the dots stand for terms of higher order in |kj_ |/«o 
and kq << G is the small parameter of the theory 
We note that t +0 > (bonding) while £_ < (anti- 
bonding). 

The quasi-particle spectrum is given by 

w - e a (k x ) -A 2 (lu + e a {k x ))- 1 - f„(k x ) = 0. (5) 

At this point we note that the RPA Green's function © 
together with JSJ bears a remarkable resemblance to the 
single electron Green's function of underdoped cuprates 
conjectured in |gj on phenomenological grounds. In both 
cases the numerator of the self energy is modified. In 
contrast, for a conventional superconductor i (kj_) would 
be expected to modify e a (k x ). 

The dispersion relationships, E a ±(k), of the quasi- 
particles are given from Eq. [SI by 

E a± (k) = t a (k ± )/2 ± VMM + t a (k ± )/2)i + A 2 . (6) 

The FS are determined by solving E a ± = 0. Doing so to 
leading order for kj^ near and G/2 yields 

[2(* x " k Fa )v F T tao? + 2t 2 aQ k ± 2 /n 2 Q = t 2 a0 - 4A 2 . (7) 

Gapless excitations, i.e. Fermi surfaces (FS), then only 
emerge when max|t a (kj_)| > 2 A. For ±t a o < 0, the FS 
pockets are electron-like, while for ±t a o > 0, they are 
hole-like. Now G a RPA = {lo + e (k))/(w - E a+ (k))(uj - 
E a _(k)). Near the FS of the pockets, e a (k x ) ~ ±A and 
E a+ (k) ~ ±t o- The effective quasi-particle weight can 
then be read off to be Z RPA ~ 1/2. Thus RPA yields 
well defined quasi-particles. 

Luttinger Sum Rule: The Green's functions and 
© satisfy the Luttinger Sum Rule (LSR) in the form 
which, though being well known theoretically, has had 
limited applications. The LSR relates the electron den- 
sity to the volume in momentum space in which G{lo = 
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FIG. 2: Electron (red and magenta) and hole (green and 
blue) pockets. The difference in size between pockets formed 
from bonding and antibonding orbitals originates from pos- 
sible difference between the corresponding tunneling ampli- 
tudes. The dashed lines represent the Luttinger surfaces 
e±(k) — 0. The thick dashed red lines are loci of the gap 
minima, e a (k) — —t a (k)/2. 



0, k) > 0. This volume is bounded by the surface where 
G(oj = 0,k) changes sign 0. The sign change can occur- 
either at an infinity of G(ui = 0, k) (Fermi surface) or a 
zero (Luttinger surface). The first possibility is denied 
for a Mott insulator. For example, the Green's function 
(0) at cj — vanishes at k Fa where taikFa) — 0, i.e. 
at momenta where the noninteracting Green's function 
had infinities. In this way the LSR is satisfied for this 
nonperturbative ladder ground state. In the presence of 
interladder tunneling, the Green's function contin- 
ues to have zeroes at kp a independent of the kj^ compo- 
nent (see Fig. 2). However, when Fermi surface pockets 
appear, the Green's function additionally changes sign 
through the newly formed infinities. Electron like pock- 
ets add and hole like pockets subtract from the total elec- 
tron density, but the LSR remains valid. This example 
demonstrates, as does the case described in 0, that in 
doped spin liquids it is generally necessary to determine 
both the Fermi and Luttinger surfaces in order to obtain 
the electron density from the LSR. It is important to 
point out that the Luttinger surface, determined by the 
zeroes of the Green's function, differs dramatically from 
the surface of minimum gap (see Fig. 2). The latter is 
often used in ARPES experiments to extrapolate to an 
underlying Fermi surface. This however leads to diffi- 
culties in the pseudogap phase of underdoped cuprates 
where the enclosed area manifestly exceeds one electron 
per unit cell, inconsistent with hole doping (e.g. see |l(1|). 
Instabilities and Doping Dependence: As was 
demonstrated in ||, the RPA solution may become un- 
stable at low temperatures. The instability is driven by 
the residual interactions between Fermi quasi-particles 
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FIG. 3: Approximation of the four-quasi-particle interaction 
by the emission of an intermediate vector boson. 



and collective modes of the spin liquid. This interaction 
may receive added strength from nesting of the Fermi sur- 
faces of particles and holes. To describe the instability 
one needs to move beyond RPA. We follow here the paper 
H and write down an effective action for quasi-particles 
interacting with collective excitations. The interaction 
comes from the diagram depicted on the l.h.s. of Fig. 3. 
This four-point function can be approximated as shown 
on the r.h.s. of this same figure, leading to the following 
effective action (p — (w,k)): 

S = \ E M-P)[" 2 - (vFk x ) 2 - A 2 ]A,(p) + 



E 

p j a 

E 

q . k 



^±(-p)G&aO#±(p) 



(8) 



r( J, -)A 3 (q)^(k)(C^U^(-k - q), 



where all fields are real and a is a charge conjugate of 
a. Aj is the bosonic field transforming according to the 
vector representation of the SO (8) group, ip± are spinor 
fields of right and left chirality, and j J are gamma ma- 
trices of the SO (8) group. In principle, there is an inter- 
action within each particle multiplet, but we neglect it 
by accepting the approximation of Fig. 3. Such interac- 
tions lead to the creation of bound states with spectral 
gap \/3A. We, however, treat these as high energy pro- 
cesses. From general considerations, supported by the 
calculation that follows, we conclude that T ~ ^/vpA. 

In coupling the ladders together, the SO (8) symme- 
try is reduced to SO(6) (provided \t+o\ = t_o| - The 
quasi-particles that transformed in an 8 dim. representa- 
tion under SO (8) now are arranged into 4 dim. spinors. 
These spinors are precisely the same that appear in the 
SO(5) theory of superconductivity ^]|. We thus expect 
the same phenomenology present in SO (5) models to be 
present here. At half-filling the coupling of the ladders 
will lead to a spontaneous breaking of the SO (6) sym- 
metry. Possible ordered states include superconductiv- 
ity (SC), antiferromagnetism (AFM), and a staggered 
flux phase (SFP). The physics of explicit SO(6) break- 
ing terms will be studied in later work. 

Moving away from half-filling introduces a nonzero 
chemical potential /i, and the SO (6) symmetry is reduced 
down to the SU(2)xU(l). The chemical potential acts 
on the vector bosons as a "magnetic field" moving the 
Cooperon down in energy. At the same time it partially 
removes the nesting such that the electron and hole pock- 
ets become unequal in size and pockets of one type may 



even disappear. In this case SC becomes the leading in- 
stability. 

We are able to provide an estimate of the supercon- 
ducting ordering temperature, T c . The RPA result re- 
mains valid at /i < A/2 so that the Cooperon has not 
yet condensed and the ground state of the single ladder 
remains unchanged. (For A sa 2/i where Cooperons do 
condense, a more sophisticated approach similar to [l2] | 
is required.) The dispersion relations are modified to 



E a (k) 



(k x - pa) 2 



E a (k + G/2) 



2mii 2m± 

(k x -PG/2) 2 



e F - [i\ 



(9) 



2 m 



2mi 



6 F -fl, 



where e F = t a0 /4 - A 2 /t aQ , P(i t i)g/4 = k Fa ± t a0 /2v F , 
m± = Ky/taQ, and my = t a o/2v' F . In two dimensions the 

density of states on the Fermi surface is p F — v ""|"' 1 = 
-jSs — . The pairing susceptibility is 

x- 1 = [(u + 2/i) 2 - (v F q x ) 2 - A 2 ] + r 2 n(^, q). (10) 

Since the interaction decays at high energies at the scale 
A, we can take it as the high energy cut-off in the po- 
larization operator. At uj,q — we have x _1 (0, 0) ~ 
-A 2 + 4/i 2 + (v F p F T' 2 /A) ln(y/(e F + h)A/T) which de- 
termines the mean field temperature of the transition to a 
superconductor with a stiffness determined by the dopant 
density: 



T c « VM*ao)+M)A 



exp 



1 - (2 M /A) 2 



const x (k v f /A) 



(11) 



This expression is valid only when t a o exceeds 2A. As fi 
increases so does T c . The T c of other possible instabilities 
(AFM and SFP) is found approximately (ignoring the 
consequences of the destruction of a (ir, tt) nesting) by 
setting /j, — in Eqn. ^2 These instabilities are thus 
exponential disfavoured. 

While the expression for T c is not valid if fi > A/2, 
we can still ask what happens to the Fermi surface for 
temperatures above any putative T c . Upon increasing 
/i > A/2, the single particle gap on the ladder decreases 
but never vanishes Q and so the ladder Greens function 
retains its zeros. It is these zeros that prevent the RPA 
electron pockets from merging together (see Fig. 2). We 
thus do not expect, within the validity of the model, a 
transition to a large Fermi surface for some \i c . 

To provide an estimate of T entering the expression 
for T c , we calculate the three-point correlation function 
in the SO (8) Gross-Neveu model using the formfactor ap- 
proach. To evaluate this correlator, we insert a resolution 
of the identity between fields, reducing the correlation 
function to a sum over matrix elements. Keeping only 
matrix elements involving single particle states we find 
upon Fourier transformation, 
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F 

(0\A a \a, Pl ) 



<o > -pi|^°|0) 
E 1 (E 1 +u>i) 



(0|V?|a, -Pi) (a, |o, -Pi) , (OICIA "Pi - P2)(/9, -Pi - p 2 |V^h -Px) 



E 2 (u)2 — E 2 ) 



E 12 (ui + lu 2 + E 12 ) 



+ 



(a, Pl \r_\p,- P2 )(p,- P2 \^\0) (a,p 1 \Tp%\0,pi+p 2 )(P,p 1 +p 2 \r_\O) 



E 2 (lo 2 + E 2 ) 



Ei 2 )(loi + lu 2 — E\ 2 ) 



(12) 



(0|^|a,p 2 )(a,P2|A a |/3,p 1 +p 2 )(/3 ! p 1 + P2 |^|0) , <0|^, -pi - p 2 )(P, -Pi - p 2 |A a |a, -p 2 )(a, -p 2 \^ |0) 



E\ 2 E 2 {lo 2 — E 2 )(uii 



E12) 



Ei 2 E 2 (oj 2 + E 2 )(uji + uj 2 + E\ 2 ) 



Each state is labeled by its isotopic index and momen- 
tum, p. Momentum and energy are parameterized in 
terms of rapidities, 9i via P i — Asmh(9i)/vF, Ei = 
Acosh(6» i ), and E 12 = \/vp{pi +p 2 ) 2 + A 2 . The matrix 
elements of the Fermi operators are given by 



(O\^\p,6)=Ae ±[ ^ i C ap e ±0 / 2 ; 



(13) 



9(0) 



exp 



B 



1/2- cosh 9 

dx sm 2 (x9/2ir) 



2cosh(cc/6) +e~ 7x/6 



iq xsinha;cosh(a;/2) 
where C is the charge conjugation matrix. A and B 



are related constants on the order of ~ yKJvp. For 
the Bose operators we have (0\A a \b,9) = (b,9\A a \0) = 

A B Sab, (p\A a \f]) = 0, where A B is 0{v~ 1/2 ). The ver- 
tex is then given in terms of the three point function 
via = 27rG- 1 (pi)G- 1 (p 2 )Gp 1 (- Pl ~ P2 )D a a(3 ( Pll p 2 ) 
where 



K (Pl,P2) 



-l7r/4 



A B A(Cr)ape- e ^/ 4 g(9 12 ) 



E\E 2 

e i^/4 



(14) 



^inh(8 2 ) 



{w\ — Ei)(u)2+E2) {uji+Ei)(lu 2 -E 2 ) 

-((P1,P2,V1,U2) -> (-Pl,Pl +p 2 ,-Wl,-^l - U 2 )), 



and the G' s are the corresponding non- interacting propa- 
gators. We also note that g(x) ~ — Be~^ x ^ 4 for \x\ >> 1. 
As we see, the vertex is a smooth function of momenta 
and frequencies changing with a characteristic scale A as 
written in (JSJ. This derivation justifies Eqs. I jllj l. 

In conclusion, we have constructed a toy model of a 
doped spin liquid. This model possesses a number of 
interesting features. Vis-a-vis ARPES measurements, it 
offers an alternative framework in which to understand 
the observed arcs in underdoped cuprates 0, 0] : such 
arcs may be Fermi pockets unresolved by ARPES due 
to disorder and limited experimental accuracy. It fur- 
ther suggests using an observed line of minimal gap will 
lead to overestimate of the number of electrons present 
in a band. Beyond implications for ARPES, the model 



(at half-filling) possesses an SO (6) symmetry and so en- 
compasses the same phenomenology as SO(5) models of 
superconductivity including a 7r-resonance at energy 2p. 
Away from half-filling superconductivity is preferred and 
the model is under sufficient control to provide an es- 
timate for the superconducting T c . Finally we point 
out above T c the model predicts the existence of a low 
lying Cooperon excitation. This excitation's existence 
and its concomitant near coherent propagation provides 
a means to understand the giant proximity effect seen in 
LSCO/LCO/LSCO thin films where the Josephson cur- 
rent was measured through LCO in its normal state . 
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